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Abstract. Dynamo action in fully convective stars is a debated issue that also questions our understanding 
of magnetic field generation in partly convective Sun-like stars. During the past few years, spectropolarimet- 
ric observations have demonstrated that fully convective objects are able to trigger strong large-scale and 
long-lived magnetic fields. We present here the first spectropolarimetric study of a sample of late M dwarfs 
(M5-M8) carried out with ESPaDOnS@CFHT. It reveals the co-existence of two distinct types of magnetism 
among stars having similar masses and rotation rates. A possible explanation for this unexpected discovery 
is the existence of two dynamo branches in this parameter regime, we discuss here the possible identification 
with the weak vs strong field bistability predicted for the geodynamo. 
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1 Introduction 



Stars: low-mass, Planets and satellites: magnetic fields, 



In cool stars, which possess a convective envelope, magnetism is thought to be constantly regenerated against 
ohmic decay by dynamo effect. For Sun-like stars the interface layer between the inner radiative zone and 
the outer convectiv e envelope is generally thought to play a major role in the dynamo processes (see e.g., 
Charbonneau 2010h . Since fully-convective stars - either main sequence stars below 0.35 M Q (i.e. with spectral 
type later than ~ MA) or young T Tauri stars - do not possess such an interface layer, generation of magnetic 
field in their interiors is often thought to rely on a non-solar-type dynam o. However, the precise mechanism and 
the properties of the resulting magnetic h ave been a debated issue (e.g., Durnev et al.lll993t Chabrier fc Kiiker 
2006; iDobler et all 120061: iBrownind 120081 ) . Two main complementary approcahes are successfully applied to 
study magnetic fields close to the fully-convective transition. On the one hand, by modelling Zeeman broadening 
of photospheric spectr al lines it is possible to assess the magnetic field averaged over the visible stellar disc (e.g., 
Reiners fc Basrill2006l) . This method is therefore able to probe magnetic fields regardless of their complexity 
but provides very little information about the field geometry. On the other hand, the Zeeman-Doppler imaging 
technique models the evolution of polarization in spectral lines during at least one rotation period in order to 
reconstruct a map of the large-scale component of the vector magnetic field on the stellar photosphere. 

Spectropolarimetric studies of a sample of M0-M4 dwarfs, conducted with ESPaDOnS and NARVAL, have 
revealed for the first time a strong change in large-scale magnetic topologies occurring close to the fully-convective 
boundary. Stars more massive than 0.5 M exhibit large-scale fields of moderate intensity featuring a signifi- 
cant toroidal co mponent and a stron gly non-axisymmetric poloidal component, with evolution happening on a 
timescale of 1 yr ([Donati et al. I l2008l D08). For those in the range 0.15-0.50 much stronger large-scale fields 
are observed, wh ich are dominated by the axial dipolar component and show only very limited evolution over 
successive years (jMorin et al.ll2008al lbl. M08a,b). Comparisons of these large-scale magnetic field measurements 
with X-ray activity indices or with mesaurements of the total magnetic field (i.e. at all spatial scales) derived 
from the analysis of Zeeman broadening of FeH molecular lines, suggest that fully-co nvective stars are muc h 
more efficient at generating large-scale magnetic field than partly-convective ones (D08. iReiners fc Basri|[2009h . 
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2 Surface magnetic fields of late M dwarfs 



A sample of 11 active M dwarfs with masses significantly below the fully-convective boundary (0 .08 < M+ < 
0.21 M(7) or spectral types M5-M8) has been observed with the ESPaDOnS spcctropolarimeter (jMorin et al 
20ld hereafter M10). Below 0.15 M , we observe two radically different categories of large-scale magnetic holds 



either a strong and steady almost dipolar held (hereafter SD, similar to stars in the range 0.15-0.5 M ); or 
a weaker multipolar, non-axisymmetric field configuration undergoing dramatic evolution on a timescale of at 
most 1 yr (hereafter WM). However the two groups of objects cannot be separated in a mass-rotation diagram, 
see Fig. [T] No object is observed to evolve from one type of magnetism to the other during the survey (some 
objects were observed for 4 years). In terms of large-scale magnetic field values, a gap exists between these 
two types of magnetism, with no object with 200 < By < 900 G in this mass range, see Fig. [5J Both stars 
hosting weak multipolar (WM) or strong dipolar (SD) fields have very strong total magnetic fields (2-4 kG). 
No systematic correlation is found between the type of large-scale magnetic topology and the total magnetic 
field Bi (see Fig. [2]). Hence, the two different types of magnetic field configurations are only detected when 
considering the large-scale component (probed by spectropolarimetry, and which represents 15-30 % of the total 
flux in the SD regime, but only a few percent in the WM regime) and not the total magnetic flux derived from 
unpolarised spectroscopy. This unexpected observation may be explained in several different ways: for instance, 
another parameter than mass and rotation period (such as age) may be relevant, two dynamo modes may be 
possible or stars may switch between two states in this mass range, etc. 
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Fig. 1. Mass-period diagram of fully-convective stars derived from spectropolarimetric data and Zeeman-Doppler Imag- 
ing (ZDI). Symbol size represents the reconstructed magnetic energy, the color ranges from blue for a purely toroidal to 
red for a purely poloidal field, and the shape depicts the degree of axisymmetry from a sharp star for non-axisymmetric 
to a regular decagon for axisymmetric. For a few stars of the sample iMorin et all (|2010h could not perform a definite ZDI 
reconstruction, in these cases only an upper limit of the rotation period is known and the magnetic flux is extrapolated, 
those objects are depicted as empty symbols. The theoretical fully-convective limit is depicted as a horizontal dashed 
line. Thin solid lines represent contours of constant Rossby number Ro=0.01 (left) and 0.1 (right), as estimated in 
IMorin et all l |2010h . 
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Fig. 2. Average large-scale magnetic fluxes of fully-convective stars derived from spectropolarimetric data and Zeeman- 
Doppler Imaging (ZDI), as a function of mass (left panel) and mass x P ro t 2 (right panel). Symbols are similar to 
those used in the mass-period diagram (see Fig. [lj . For stars in the WM regime symbols corresponding to different 
epochs for a given star are connected by a vertical grey line. The yellow region represents the domain where bistability 
is observed and the orange one separates the two types of magnetic fields identified (see text). 



3 Weak and strong field dynamos: from the Earth to the stars 

In this section we briefly detail one of the hypothesis that could explain the observation of two groups of late 
M dwarfs with very different magnetic properties: the fact that two different dynamo modes could genuinely 
operate in stars having very similar mass and rotation. We focus here on the weak vs strong field dynamo 
bistability, initially proposed for the geodynamo. The underlying idea is that both, magnetic fields and rotation, 
taken separately tend to inhibit convection, but that if both effects are combined the impeding influences of 
the Lorentz and of the Coriolis forces may be reduced, allowi ng convect ion to se t on at lower Rayleigh number 
and to develop on larger length scales ()Chandrasekharj[l96ll) . This led Roberts! ()l978l ) to conjecture that in a 



rapidly rotating system, for magnetic fields stronger than a threshold value, the Lorentz force would enhance 
convection and hence dynamo action, resulting in a runaway growth of the magnetic field. The corresponding 
bifurcation diagram is depicted on Fig. [3] On the weak-field branch the Lorentz force is balanced by viscous or 
inertial terms in the momentum equation, this force balance requires small-spatial scales. On the strong field 
branch, however, the magnetic field strength is set by a balance between Lorentz and Coriolis forces, which 
requires larger spatial scales, this is the magnetostrophic regime. A similar bifurcation diagram, but based on 
the fact that ma gnetic buoyancy would be negligible close to the dynamo onset has been proposed for stars by 
Weiss fc Tobial (l2000l) . 



We now briefly discuss the id entification between WM (SD) magnetism and weak-field (strong-field) dynamo 



regime, the reader is referred to iMorin et al. (2011) for a more detailed discussion. First, in the SD regime the 



magnetic field strength has to be compatible with a Lorentz-Coriolis force balance. We note that this balance 
is valid spatial large-scales for which the Coriolis term is predominant over inertial terms in the momentum 
equation, in qualitative agreement with the observation that only the large-scale component of the magnetic 
field is exhibits a bimodal distribution. This magnetostrophic force balance roughly corresponds to an Elsasser 
number of order unity, i.e. : 

A = -^-~l, (3.1) 
Plirjil 

where B is the magnetic field strength, p the mass density, p the magnetic permeabilit y, rj the magnetic diffusivity 
and Q the rotation rate. With a few assumptions described in Morin et al. ( 201 lh . we find that the order of 



magnitude of the expected magnetic field strength on the strong field branch is set by: 



B. 



of 



m m m m m) * 

•ill „„2 , 



Where ry© is a reference value for the magnetic diffusivity in the solar convection zone, and 77 rB f = 10 cm s 



Takin g stellar radius and luminosity for the stellar mass in the range 0.08 — 0.35 M Q from Chabricr & Baraffe 
( 1997 ) main sequence models, and tjq in the range 10 11 — 3 x 10 12 cm 2 s _1 (e.g., Riidiger et al.l 201 ll) . we 



240 



SF2A 2011 



derive surface values in the strong field regime in the range 2-50 kG, compatible with the order of magnitude 
of measured By values. More conclusively, the gap in terms of By between the two branches depends on the 
ratio of inertia to Coriolis force in the momentum equation and can estimated with: 

= Ro 1 ' 2 , (3.3) 

B sf 

which is of the order of 10 _1 for stars of our sample in the bistable domain, in good agreement with the typical 
ratio of large-scale magnetic fields measured between the SD and WM groups of stars (see Fig. [2]). 




Fig. 3. Left: Total magnetic fluxes of fully-convective stars measured from unpolarised spectra of FeH lines. The 
values are taken from iReiners et al.l (|2009T ) and iReiners fe Basril (|2010h . whenever 2MASS near infrared luminosities 
l|Cutri et alj|2003 l and Hipparcos parallaxes (?) are available to compute the stellar mass from the lDelfosse et all (2000) 
mass-luminosity relation. Whenever spectropolarimetric data are available the properties of the magnetic topology are 
represented as symbols described in Fig. [2] Else small blue symbols are used, upward (downward) triangles represent 
lower (upper) limits. Stars with Bf < 1 kG are not in the saturated regime of the dynamo (slow rotators) and are 
shown for comparison. Right: Anticipated bifurcation diagram for the geodynamo (adapted from lRobertsl[l988l b The 
magnetic field amplitude is plotted against the Rayleigh number. The bifurcation sequence is characterised by two 
branches, referred to as weak and strong field branches. The yellow and orange regions have the same meaning as in 
Fig. [2] Ra c is the critical Rayleigh number for the onset of non-magnetic convection. The weak field regime sets in at 
Ra m , and the turning point associated with the runaway growth corresponds to Ra = Ra r . 



4 Conclusions 

Brief summary; hysteresis; impact on stellar formation/evolution 
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